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Abstract 

We show that photon coincidence spectroscopy can provide an un- 
ambiguous signature of two atoms simultaneously interacting with 
a quantised cavity field mode. We also show that the single-atom 
Jaynes-Cummings model can be probed effectively via photon coinci- 
dence spectroscopy, even with deleterious contributions to the signal 
from two-atom events. In addition, we have explicitly solved the eigen- 
vectors and eigenvalues of two two-level atoms coupled to a quantised 
cavity mode for differing coupling strengths. 

1 Introduction 

Photon coincidence spectroscopy (PCS) has been introduced as a tool 
for probing the spectrum of the Jaynes-Cummings ( JC) model corre- 
sponding to a single two-level atom (2LA) coupled to one quantised 
field mode of an optical resonator. Such a system is realised via atomic 
beam cavity quantum electrodynamics (CQED) P HI H U E 13 El El 
110] . In two-photon coincidence spectroscopy (2PCS), a sparse atomic 
beam traversing the optical cavity is driven off-axis by a bichromatic 
driving field, and two photons in the output field of the cavity are 
deemed to be coincident if they arrive at the photodetectors within 
some short time interval [01 EJ. During the 2LA's passage through 
the cavity, the electric dipole coupling strength g(r) varies with its 
position r, and the coupling strength peaks at an antinode along the 
cavity's longitudinal axis. The 2LAs are assumed to be sufficiently 
slowly moving so that the system can be treated as a JC stationary 
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system with a coupling strength g varying according to the proba- 
bility distribution P(g). Thus, atomic motion is responsible for an 
effective inhomogeneous broadening of the spectrum for the coupled 
atom-cavity system. However, the bichromatic driving field selectively 
excites a subensemble of coupled atom-cavity systems (with particular 
values of g) to the second couplet of the JC ladder. The JC spectrum 
is then probed via two-photon coincidence spectroscopy by measuring 
the rate of two-photon coincidence events, which is referred to as the 
two-photon coincidence rate (2PCR) 

A bichromatic driving field has been used in optical CQED 4 , 
but photon coincidence measurements are yet to be implemented in 
experiments. More recently experiments have succeeded in controlling, 
to some extent, single-atom motion in the cavity but complete 
confinement of a single atom in a cavity [12111] with P(g) oc 5(g-g Q ), 
i.e. a 'frozen' atom, is yet to be achieved. 2PCS is designed specifically 
to overcome inhomogeneous broadening due to atomic motion, thereby 
permitting the nonlinear component of the JC spectrum (specifically 
the second couplet of the JC ladder and its two photon de-excitation 
events) to be optically probed. 

2PCS has been shown to be effective for overcoming the challenges 
of inhomogeneous broadening due to atomic motion. However, a sec- 
ond challenge is atom number fluctuations. As the atoms arrive at 
the cavity randomly in an atomic beam, the number of atoms in the 
cavity changes: most of the time no atom is in the cavity; some of the 
time there is one atom; less frequently there are two or more atoms in 
the cavity. 

In quantum trajectory simulations of 2PCS, multi-atom effects 
have been included in the simulations [5], and the nonlinear regime of 
the JC spectrum may be observed even in the presence of multi-atom 
contributions to the photon statistics. Of course a sparse beam of 
2LAs is required to ensure that multi-atom effects do not overwhelm 
the photon statistics. Rigorous bounds on the mean atom number 
in CQED have been established [Tl|- Here we investigate in detail 
the contributions of multi-atom effects to the extraction of an un- 
ambiguous signature of the JC spectrum. We show additional peak 
structure that arises in the two-photon coincidence spectrum due to 
the two-2LA contributions, but we also show that these effects can 
be negligible for a sparse beam. An effect of multi-atom contribu- 
tions to the 2PCR is to raise the background 2PCR, thereby making 
2PCS only marginally more difficult than if two-atom effects could be 
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avoided altogether. 

Furthermore, we consider the system of two 2LAs coupled to a 
single cavity field mode. For the sparse atomic beam passing through 
the cavity, we show that an off-axis bichromatic driving field can 
be used to observe the spectrum of the coupled system consisting 
of two 2LAs and one quantised cavity field mode. By choosing the 
two frequencies of the bichromatic driving field judiciously, one may 
extract the spectrum for various choices of coupling strengths g\ and (72 
for atoms 1 and 2, respectively. The case g\ = 52 corresponds to 
probing the two-atom version of the Tavis-Cummings model |15| . 

We develop the Hamiltonian for multi-atom systems in sectional 
In section E3 we develop 2PCS as a tool for probing the spectrum of 
two atoms coupled to a single quantised cavity field mode, and we 
conclude in section @J 

2 Mathematical background 

In the typical CQED experiment, which tests for quantum field effects, 
an atomic beam is directed through a cavity [U El and is excited by 
a driving field during the passage of the atom. In contrast to mi- 
crowave cavity quantum electrodynamics El) where the atomic 
lifetime is large compared to the passage time through the cavity, the 
lifetime in an optical system is quite short, hence the need for opti- 
cally driving the system during the passage of the atom through the 
cavity. Whereas current experiments employ an on-axis driving field 
(the laser field is directed into the cavity mode through one mirror), 
2PCS is dramatically improved for an off-axis driving field (the atoms 
are driven directly during passage through the cavity by a bichromatic 
driving field), and photon coincidences are detected in the cavity out- 
put field. 

2.1 The Hamiltonian 

The Hamiltonian for the cavity field is given by 

H<® = htoa^a , (1) 

for n = a^a the photon number operator, and u is the angular fre- 
quency of the cavity mode (assumed to be equal to the transition 
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frequency for the 2LA). If there are N 2LAs in the cavity, the Hamil- 
tonian for the m th 2LA is generalised from iJw of equation Q to 
yield 

H {m \g m ) = huai^+ih^a^-g^aa^), (2) 

with g m = g(r m ) the coupling strength between the m th atom and 
the cavity field. Although the atom is moving, we assume that the 
motion is slow so that time-dependence of the Hamiltonian can be 
ignored. Atomic motion is then accounted for by allowing g m to be 
a random variable [HUHj- Taking g m = \g m \ exp(i6 m ), and applying a 
unitary transformation U^ m \g m ) = exp(i6 m ai m ' ') to equation (j2}, the 
Hamiltonian simplifies to 

[U^]\g m )H^\g m )U^\g m ) = H^\\g m \). (3) 

Thus, we use the simpler Hamiltonian (J3J) for treating the m th atom 
and assume that g m is real and positive. The 2PCR is a function of 
the two-photon correlation 

w^\t) = (:n(t )n(t + t):) (4) 

for an unimportant 'initial' time to, and '::' imposes normal ordering. 
The time interval for coincidences is assumed to be short so that the 
2PCR is approximately proportional to [HI El 

u> (2) = (: n 2 :). (5) 

The atomic beam is slowly moving, and the time dependence of g m 
can be ignored by instead considering a distribution P(g m ) for the 
varying coupling strength of the m th 2LA. For the multi-atom sys- 
tem, each atom experiences a varying coupling strength. We assume 
a cut-off of N atoms in the system; more than N atoms introduces a 
negligible effect to measurable quantities. We can express the coupling 
strength for the N atoms by the multivariable coupling-strength dis- 
tribution P{g) for g = (gi,g%, . ■ ■ ,gjsr)- The Hamiltonian for N 2LAs 
in the cavity is given by 

N 

H N (g) = H^ + Y,^ m \g m ). (6) 

m=l 

Under the unitary transformation 

U N (g) = UW( 9N ) ® U^ N - l \g N ^) ■ ■ ■ U (2 \g 2 ) £/ (1) (<7i), (7) 
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equation © reduces to 

N 

U ] N (g)H N (g)U N (g) = Wo + £ hue™ + ihg m (at*™ - aaf } ) , 

m=l 

(8) 

where we assume that g m is always real and positive. The assump- 
tion is valid because (t) in equation is invariant under unitary 
transformations by Ujy(g). 

The case Hm<n {.(})■, corresponding to fewer than N atoms in the 
cavity, is implicit by setting one or more values of g m to be zero. For 
N = 1 the multiatom Hamiltonian © reduces to the JC Hamilto- 
nian ^Bl with 'dressed states' 

|0> = |0) cav |g> (9) 

and 

i 

\n)± = -j= (\n - l) cav <g> |e) ± i|n) cav ® |g)) . (10) 

Here {|g), |e)} are the ground and excited state of the 2LA. The N = 2 
case is developed in detail in Appendix A. Instead of a ladder consist- 
ing of a ground state as a singlet and higher-order states as doublets, 
there is a singlet, then a triplet (for n = 1) and then quadruplets (for 
n > 1). This assumes that the two atoms experience different coupling 
strengths gi ^ g2- For g\ =52, the two-2LA Tavis-Cummings model 
is obtained [To] . 

2.2 The master equation 

The Hamiltonian evolution (JSJ does not account for either irreversible 
dynamics or contributions from the bichromatic driving field £{t) = 
Si exp (— iuit) + £2 exp (— wit) (with the resonance condition oj\ = 
oj — gf). To include all these effects, we construct the quantum master 
equation. The Born-Markov approximation is applied to both radia- 
tion reservoirs: the reservoir for the field leaving the cavity and the 
reservoir for direct fluorescence of the 2LA from the sides of the cavity. 
The cavity damping rate is k, and the emission rate into free space 
is j m , where j m is the inhibited spontaneous emission rate for the 
m th atom. The master equation [7] can be expressed as p = Cp for C 
the Liouvillean superoperator and p the density matrix for the system. 
More specifically the Liouvillean superoperator can be expressed as the 
sum of a time-independent Liouvillean operator, a time-dependent 
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Liouvillean operator and a 'jump' term. In the rotating frame, the 
time-independent non-Hermitian effective Hamiltonian is 



N 



(m) , 



(11) 



m=l 



with 



H 



(m) 
eff 



( (oj — uj\) a'a — iKa'a, 



i \ ( m ) i ■ t t ( m ) 



(m)\ 



if m = 



m > 0. 



(12) 



-TH^)-i( 7m /2)4 m) a W 
The monochromatic driving term is given by 

TH(£i)=tfi(4 m) -ffl m) ). 
Thus, the Liouville operator for the iV-atom cavity system is 

AT 



(13) 



C P (g) 



n 



m=l 



£ - f + 2>fc c/)p + (14) 



with time-dependent driving term 



T>(t;6)p 



T(£ 2 e- M ),p(g) 



and jump term 



N 



Jp = 2Kap(g)a) + £ p(<f)<r 



(15) 



(16) 



m=l 



We assume that 7 = j m is identical for all atoms. For atom number 
fluctuations, the 2PCR is the sum of 2PCRs for each case: no atom, 
one atom, two atoms, and so on. We calculate the 2PCR for each 
case and add each 2PCR component to obtain the resultant w^ 2 \ 
the 2PCR of equation ®, up to an unimportant scaling coefficient. 
However, the contribution to the 2PCR must be weighted according 
to the probabilities {p m }, with p m the probability for m atoms being 



in the cavity. Thus, the total 2PCR is = Em=o^ 



,w 



(2) 
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2.3 The density matrix 



The master equation is used to solve the density matrix p in the 
long-time limit, well after transients have disappeared. However, the 
bichromatic driving field results in an oscillatory solution for p; a time- 
independent stationary solution is not obtained for this case. The 
long-time limit for p can be solved by employing a Bloch expansion in 
terms of the bichromatic driving field frequency difference 5 = uj% — uj\ , 



with pk(t, g) time-dependent matrices. In the long-time limit, p\ « 0, 
and Pk(t,g) — * Pk(g) is time-independent. As the photocount integra- 
tion time is expected to be long compared to the frequency 5, it is rea- 
sonable to approximate p(t,g) by truncating expansion (|17|). For pho- 
tocount integration times larger than 5~ 1 , we truncate (|17|) at k = 0, 
yielding p(t,g) ~ Po(g), which is time-independent. Approximating 
the time-dependent density matrix by the first (time-independent) 
term in the Bloch expansion is consistent with single-atom 2PCS stud- 
ied in reference [7]. 

The resultant density matrix for the A-atom cavity system is given 



where g ma .x is the coupling strength at a cavity antinode, Fg max is the 
effective lower bound cut-off for the coupling (0 < F < 1) and P(g) = 
P(gi)P(g2) ■ ■ ■ P{gN) is the iV-atom coupling strength distribution, 
with two typical plots of P(g) for one atom depicted in figure Q for k 
the cavity loss rate [Hl[7j. The effect of averaging over P(g) is an inho- 
mogeneous spectral broadening, and we use the overbar ~ notation to 
denote averaging over P(g) to account for inhomogeneous broadening. 
This broadening is due to atomic position variability. One distribution 
(dotted line) corresponds to the case of a uniformly distributed atomic 
beam entering the cavity and the other (solid line) to an atomic beam 
initially passing through a rectangular mask [2j . In both cases we have 
assumed a single-mode cavity supporting a TEMoo mode and atomic 
motion perpendicular to the longitudinal axis of the cavity. 
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(17) 



k=0 



by 




(18) 



7 



3 Photon— coincidence spectroscopy for 
two atoms in the cavity 

Although PCS has been developed to probe the nonlinear portion of 
the (one-atom) JC spectrum (2PCS probes the second couplet and 
multiphoton coincidence spectroscopy [Hj probes higher levels), in this 
section we show that 2PCS can also probe the nonlinear portion of the 
spectrum for two-atom events. Interpreting the 2PCR for the two- 
atom case is somewhat more complicated than for the JC spectrum 
due to the presence of inhomogeneous broadening for both atoms; 
methods for interpreting the 2PCR peaks must therefore be more so- 
phisticated. In this section, we describe the method for performing 
PCS to extract a quantum field signature for two atoms in the cavity. 

3.1 Excitation of two atoms coupled to a sin- 
gle cavity mode 

Two-atom dressed states are discussed in detail in Appendix A, and 
this description helps to understand the efficacy of PCS for two-atom 
CQED. The strong-coupling regime implies that considering atoms 
and the cavity mode separately is not helpful: dressed states are the 
preferred description. 

The linear portion of the Hamiltonian spectrum consists of the 
ground state and the one-quantum triplet. This triplet reduces to a 
doublet for the case that both atoms experience the same coupling 
strength g\ = g2 15 and is the system being probed in some normal- 
mode (or vacuum Rabi) splitting experiments which involve multiple 
atoms in the cavity 00 E|. Such experiments have been described 
well by semiclassical theories [TJ2|. The purpose of 2PCS is to extend 
beyond the linear regime to where quantum field theory is essential for 
describing observations. Thus, for the two-atom case, 2PCS should 
be designed to probe the first quadruplet of states (or higher, as an 
analogy to multiphoton coincidence spectroscopy [S]). 

The singlet, triplet and first quadruplet of dressed states for the 
system consisting of two 2LAs coupled to a single cavity mode are 
depicted in figure b) . Two-photon excitation to the first quadru- 
plet is developed by analogy with two-photon excitation of the JC 
system to the second doublet, depicted in figure Efa), and discussed 
in references E] . 
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For the (single-atom) JC system in figure Ufa), we depict the 
two-photon excitation scheme from the ground state |0) to the first 
two excited couplets \n) e (n £ {1,2}, e G {—>+}) of the dressed 
states. This single-atom system is a useful guide to understanding the 
two-atom scheme, and we briefly explain the single-atom case here. 
The challenge is to overcome inhomogeneous broadening of the cou- 
plets |1) £ and |2) £ is 2% max and 2\^2hg mSLX , respectively, due largely 
to fluctuations in the atomic position. Two two-photon excitations 
to the second couplet of the (one-atom) JC system are depicted for 
a bichromatic driving field with one component of amplitude £\ and 
frequency uj\ and the other with amplitude £2 and frequency 102- The 
excitation pathway on the right of figure [2^a) {uj\ then 0J2) excites 
resonantly from |0) to |1)_ and then may excite resonantly from |1)_ 
to either of the states |2)±. The excitation pathway on the left {lo2 
then uji) excites resonantly from |0) to |1) + to |2)_ for g = (\/2 — 

For 2PCS, the signature of entanglement is obscured by the back- 
ground 2PCR due to two 012 photons contributing to excitation to 
the second couplet via nonnegligible, off-resonant transitions. The 
method for overcoming this problem is called 'background subtrac- 
tion'. In this technique, the experiment is repeated twice, once with 
the bichromatic driving field and again with the fixed field turned off. 
The difference between these two 2PCRs shows all the features of the 
desired 2PCR without the deleterious effects of two 0J2 transitions to 
the second couplet states [HI Ej • 

Now we consider the case of two 2LAs atoms coupled to a single 
quantised field mode, which is depicted in figure Efb). Two-photon 
excitation occurs from the ground state |0) via the triplet {1)^ to 
the first quadruplet |2) ££ , (with t] G { — ,0,+} and {e,e'} £ {—,+}) 
with gi 7^ g2 assumed (the case g\ = 52 is sufficiently unlikely to be 
ignored). Two two-photon excitation pathways to the second couplet 
are depicted for a bichromatic driving field with one component of 
amplitude £\ and the other with amplitude £2- The excitation path- 
way on the right (wi then 002) excites resonantly from |0) to |1)_ and 
hence may excite resonantly to any state of the set {|2) ee /}. The ex- 
citation pathway on the left (uj2 then u)\) excites resonantly from |0) 
to 1 1 ) ^ to |2) ££ /. By analogy with one-atom 2PCR, we can employ a 
spectral hole-burning approach for excitation to a subensemble of the 
first quadruplet and thereby operate in the nonlinear regime. 
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3.2 Two— photon coincidence spectrum 

For 2PCR of the JC system, the bichromatic driving field is used to 
excite selectively to subensembles of the inhomogeneous broadened 
first couplet of the JC ladder, as depicted in figure E^a). The second 
photon excites to a state of the second couplet only when the resonance 
condition for u>i+u> 2 equals 2u±\/2g for the selected subensemble with 
coupling strength g; this resonance signature provides unambiguous 
evidence of a quantum field effect through the "v2" ■ 

For 2PCR of the two-atom plus cavity mode system, selecting 
a particular subensemble corresponding to a fixed g = ((71,92) is not 
possible. To see how this complication arises, consider the case that ui\ 
is tuned such that direct excitation occurs from |0) to which is 
given by equation (j22j), for certain choices of g\ and g 2 . The energy of 
state |1)_ corresponds to the eigenvalue A_ of equation (|2*H|) . Choosing 
a particular oj\ fixes g = y/ g\ + g\ according to the expression A~ = 
ijj — g but does not fix g\ and 92 separately. Therefore, choosing a 
particular uj\ does not select a unique subensemble but rather a class 
of subensembles corresponding to the constraint g = to — X~ = to — uj\. 

This large class of subensembles is not too difficult to manage, 
however. We are particularly interested in factorisable distributions 
of the type P(g) = P(g 1 ,g 2 ) = P(gi)P(g 2 ); i.e. P{g m ) is identical 
for each independent atom, and P(gi) appears as shown in figure ^ 
Choosing a fixed u>i constrains the selected subensembles to values 
of g in the domain of interest, namely large g, and effectively reduces 
contributions such as high g\ and very low g 2 that would negate the 
desired 2PCR peaks. 

For 2PCR of the two-atom-plus-cavity system, the 2PCR is en- 
hanced, not for specific g\ and g 2 but rather for a range of g± and g 2 
constrained only by g being determined by the choice of uj\ (and we un- 
fortunately cannot fix both g\ and g 2 in experimental circumstances). 
This means that the expected 2PCR peaks are broadened by this lim- 
ited control of g\ and g 2 . The 2PCR for two atoms in the cavity, with 
strong coupling to the mode and a weak bichromatic driving field, is 
plotted in figure 01 

As an example, we consider the outlying strong peak at 5 = 2.414 
(peak viii of figure Efb)). As we will show, this peak is due to the 
excitation pathway |0) < — ► |1)_ < — ► |2) ++ . The frequency differ- 
ence for |0) < — > |1)_ is A - , as discussed above, and the transition 
for |1)_ < — > |2)++ is A^ + — A - . The expression for normalised de- 
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(19) 



u — U>i 



could be plotted as a line in the three-dimensional space spanned by 
the gi, 52 and 5 axes. We do not present such a plot here but rather 
note that the trajectory intersects the observed peak value 5 = 2.414 
(peak viii of figure El) for g\ = <? max arid 52 = or vice versa. This 
intersection tells us that the peak at 5 = 2.414 is really just a single- 
atom JC system peak because the second atom is effectively decoupled 
(negligible 52)- Setting the frequency u>\ to X^ 1 is therefore not the best 
way to probe two-atom effects because X^ 1 depends on g but is not 
dependent on g\ and g<i separately. This has the effect of producing a 
peak for the case that 1 1)± reduces (for g\ = g,g2 = 0) to 2~ 1 / 2 (|0)|e)± 
|l)|g))|g). Atom 2 remains in the ground state. Similarly, setting u>i 
to A has the problem that it is completely independent of uj\ and u>2- 
The dilemma of not having a good choice for u\ to excite from |0) 
to |1)-|- or I l)o is solved by setting u>2 to X^ or A . However, u>2 is 
the scanning frequency. The proposal here, though, is first to as- 
sume that ll>2 is set to the frequency A - , which excites |0) < — > |1)-- 
Then uj\ = X ^ — A - , which is appropriate for the transition |1)_ < — > 

1 2) (.. This counter-intuitive choice (letting the scanning frequency 

excite the first transition rather than the second) enables peak i of fig- 
ure El to be understood. The 2PCR would then be enhanced for 0J2 = 
A~ with uj\ fixed. The peak corresponding to this excitation pathway 
(depicted in figure EJ is observed at 5 = —1.345 (peak i of figure EJ). 
Substituting this value of 5 into the expression for 5 (g\, 52) obtained 
in Appendix A yields the two solutions gi/n = 64.9 and 52/^ = 45.5 
and the reverse. The peak at 5 = —1.345 has its most significant 
contribution from the case that both atoms contribute to the 2PCR 
signal. 

This analysis helps to understand the 2PCR peak structure de- 
picted in figure El but a detailed analysis is afforded by the method of 
suppressed transitions developed in reference ^U]. In this method we 
identify the specific transitions contributing to each peak. 

3.3 Method of suppressed transitions 

We apply the method of suppressed transitions to understand figure El 
and to identify the peaks corresponding to genuine quantum field ef- 
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fects for two atoms in the cavity. In the method of suppressed tran- 
sitions ^U], we modify the Liouvillean superoperator C in numerical 
simulations by artificially eliminating particular driving terms respon- 
sible for certain peaks. The isolation of specific transitions is obtained 
by the following procedure. For example, let us consider the influ- 
ence on the 2PCR of the |0) < — > transition. We can write the 
effective Hamiltonian (jllj) for two 2LAs as a matrix in the dressed- 
state basis. To isolate the influence of the |0) < — ► |1)_ transition, 
we can set the matrix element (0|Y(£i)|l)_ and its complex conju- 
gate _(l|Y(£i)|0) to be zero where Y(£i) is the (Hermitian) driving 
term T(£i) = TZ=l T (m) (£i) for ^ (m) {£i) defined by equation (Ety. 
In addition the matrix element (0|Y(£2 exp (— z<5i))|l)_ and its com- 
plex conjugate can both be set to zero. The jump term J is not 
modified because only the driving terms and their effects are of con- 
cern in this analysis of suppressed transitions. 

In subsection 13.21 we have shown that transitions from |0) to |1}„ 
{rj G {—,0,+}) via an u± photon is not the best choice for probing 
the two-2LA cavity system as these excitations lead to single-atom 
JC system peaks. The artificial removal of these transitions by the 
method of suppressed transitions eliminates the (single-atom) JC sys- 
tem peaks in the simulations and thereby helps to identify those peaks 
that are primarily due to two-atom cavity events. The method of sup- 
pressed transitions is a mathematical tool for interpreting peaks in the 
two-photon spectral structure, not a physical process. 

We begin by setting both the matrix element (0|Y(£i)|l}_ and 
its complex conjugate to zero as this transition is present in every 
single-atom peak, and suppressing this transition is an excellent first 
step to identify peaks associated with two atoms coupled to the cav- 
ity mode. The result is depicted in figure Efa) without and Efb) 
with background subtraction. We observe that eliminating these two 
matrix elements causes a dramatic reduction of peaks ii, Hi, vi, vii, 
viii with the labelled peak structure of figure Gth) replicated in fig- 
ure [5f b) as a dotted line. The reason for these reductions is that 
the excitation paths responsible for these peaks all have a signif- 
icant contribution from the |0) < — ► |1)_ transition, employing a 
photon of frequency uj\. Suppressing (0|Y(£i)|l)_ and its conjugate 
only supplies a partial interpretation of the peak structure. A further 
understanding of these peaks is obtained by imposing the condition 
that Y(£2 exp(z<5i))|2) ++ and its complex conjugate are zero. As 
shown in in figure EJ photon U2 dominates the |1)_ < > |2) ++ tran- 
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sition responsible for peak viii. Thus, the wi-driven excitation path- 
way |0) < — > followed by an u^-driven transition to the |2) ++ 
state is responsible for peak viii as shown in figures EJb) andHJb), re- 
spectively. This method of suppressed transitions validates the analy- 
sis of peak viii in subsection l3.21 where the excitation pathway |0) < — > 
|1)_ < — ► 1 2)++ was suggested as being responsible for the existence 
of this 2PCR peak. 

As another example, we consider the transition |1)_ < — ► |2)_| 

By eliminating this transition, peak vii (and vicinity) are reduced as 

shown in figure E[d). Thus, the |0) < — > |1)_ < — > |2)_| pathway, 

induced by an oj\ photon followed by an IV2 photon, is partially re- 
sponsible for the increase of the 2PCR at vii (and vicinity). Figures El 
and|B[f) show clearly that an wi-driven |0) < — ► transition and 
an (^"driven |1)_ < > |2) transition are responsible for the back- 
ground of peaks iv, v and vi but not for the actual peaks themselves. 

Setting _(1|T(^2 exp (i5t))\2) |_ and its complex conjugate to zero 

eliminates peaks ii and Hi as shown in figures HJg, h). Clearly, peaks 
ii and Hi are primarily due to a resonant |0) < — ► |1)_ transition via 

an oj\ photon, followed by an excitation via an L02 photon to |2) y. 

state. Peaks ii and Hi are the dominant peaks because, apart from 
the fixed field with frequency u\ driving |0) < — ► |1}_ , the scanning 
field also drives |0) < — ► |1)_ via an off-resonance transition; thus, 
higher populations of the |2) |_ is expected. 

Thus far, we have identified those peaks that are due to single- 
atom events. The remaining peaks, particularly iv, v and vii, are the 
ones of interest because they are primarily due to the two-atom cav- 
ity signature, and this is because the transition |0) < — * (1)^ state is via 
ana>2 photon. To explain peaks iv, wand vii, we set (0|X(f 2 exp(—i5t))\l) 
and r] {l\T(£i)\2) ee i, and their complex conjugates, all to zero. For ex- 
ample, setting (0|T(^2 exp(— i6t))\l}+ and its complex conjugate to 
zero eliminates the peaks at v, vi and vii as shown in figure Efb). 
With this result and figure Etb), we observe that the peak at vii 
is due to an cj2 - driven |0) < — > |1)+ transition followed by an oj\— 

driven |1) + < — > |2)_| transition. Fixing +(l|Y(£i)|2) and its 

complex conjugate to zero reduces the peak at vi as shown in fig- 
ure Efd). Thus, the small bump, namely peak vi (to the right of 
peak v), is due to an ^-driven excitation |0) < — > |1)+ followed by 

an wi-driven |1) + < — > |2) transition. As shown in figure Etb) and 

figure EJf), the removal of the u^driven transition |0) < — ► |1)+ and 
the wi-driven |1) + < — > |2) |_ transition yields a diminished peak 
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v. Peak iv is diminished if the transitions (0|T(£2 exp(— i5t))\l)o, 
o(l|T(£ 2 e3cp(-i<5t))|2)__ and (1|T(£ 2 exp(— *5t))|2>_ + , together with 
their complex conjugates, are suppressed by being set to zero as shown 
in figures [7(d), [H(b) and IHJd). Clearly, peak to is due to an uj2~ 
driven resonant pathway |0) < — > |l)o, followed by wi-driven transi- 
tions | l)o < — ► 1 2) |_ and |l)o < — ► |2) Finally, peak z vanishes if 

we eliminate the W2-driven resonant pathway |0) < — * |1)_ followed by 
the wi-driven |1)_ < — ► |2) |_ transition (figures [7_(e, f) and[B(e, f)). 

The clearest unambiguous spectroscopic feature of the two 2LA- 
cavity system, which cannot be accounted for by the (single-atom) 
JC-spectrum, is peak i. This peak is almost visible in figure EI a) and 
clearly visible after background subtraction as shown in figure 13(b). 
In addition, the peak is much more pronounced than the other rele- 
vant 2PCR peaks at iv, v and vii. Thus, this peak provides an ideal 
signature of quantum field effects in two-atom cavity systems. 

Furthermore, the method of suppressed transitions allows us to 
identify the deleterious two-photon pathways that destroy the quan- 
tum signature of single-atom and two-atom peaks (single-atom peaks 
ii and Hi and two-atom peaks i, iv, v and vii). These pathways are 
responsible for the difference between figure 0(a) and 0(b). Back- 
ground subtraction eliminates precisely those deleterious pathways. 
Figures E(g) and [7(e) clearly show that peaks %, ii and in are strongly 
affected by the |0) < — > |1)_ < — ► |2) |_ pathway via two u>2 pho- 
tons. The second major undesired pathway washes peaks iv, v, vi 
and vii away. This pathway is primarily due to the W2-driven transi- 
tion |0) < — ► 1 1)+ (as shown in figure [7(a)), followed by an ^-driven 
transition |1) + < * |2) ++ (not shown). 

4 Conclusion 

For an atomic beam the number of atoms in the cavity mode is a ran- 
domly varying quantity. For a sparse beam, most of the time there 
is effectively no atom in the cavity. A large portion of time there is 
effectively one atom, although, of course, its position is a randomly 
varying quantity. There is also a contribution to the photon coinci- 
dence spectrum due to two or more atoms simultaneous interacting 
with the cavity mode. In section |31 we have assumed precisely two 
atoms in the cavity but left the positions random. In this section 
we include the signal for one and no atom events in the cavity with 



14 



appropriate weightings for a sparse beam. We show that the 2PCR 
reveals genuine quantum field effects for two-atom events. 

We let p n represent the density matrix for n atoms in the cavity, 
averaged over the atomic positions. This matrix p n is a generalisation 
of equation (|18|) where the subscript refers to the number of atoms in 
the cavity. For the sparse atomic beam, we will ignore p n for n > 2. 
The case n = is not significant because off-axis driving ensures that 
photon coincidences do not arise when there is no atom present. Here 
we make the reasonable assumption that P1/P2 = 9. That is, the 
probability of having one atom in the cavity is nine times more likely 
than having two atoms. This is compatible with the sparse beam 
assumption. The result for p\jp2 = 9 is shown in figure HOT a. b). 

We consider two high value of <?/ (= 63k) as discussed 

in subsection 13.31 (<?/ = 63k was assumed in a study of multipho- 
ton coincidence spectroscopy for the JC model [S]), and a low value 
of gf (= 9k), which has been the subject of other studies of photon 
coincidence spectroscopy [HI El El- The case of a high value of <?/ cor- 
responds to strong coupling and is depicted in figure ITUT a): the low 
value case is depicted in figure ITr¥ b). 

For gf = 63k, we observe that peak i, which corresponds to a 
genuine quantum field effect with two strongly coupled atoms in the 
cavity, is still quite pronounced, despite the strong signal from single- 
atom effects. Even for P1/P2 = 9, the two-atom events contribute a 
strong difference-2PCR signal. Therefore, we can assert that two- 
photon coincidence spectroscopy is an effective means for extracting 
the quantum field signature of two atoms in a quantised field cavity 
for a sparse atom beam. 

On the other hand, we explore the case of relatively low coupling, 
namely gf = 9k. This case is depicted in figure HUf b) and shows 
that the single-atom 2PCR spectrum (depicted as a dotted line) is 
not modified in any substantial way by the two-atom contributions. 
Therefore, the employment of 2PCS as a probe of quantum field effects 
in the (single-atom) JC model is sound for a sparse atomic beam. Of 
course it was assumed in references (HI E] that multi-atom effects were 
not significant factors in modifying the ideal single-atom two-photon 
spectrum. However, here figure ITf¥ b) reveals quite clearly how much 
of an effect arises from two-atom events and how small this effect is. 

On the one hand, for strong coupling, the two-atom effect is large, 
and peak i, in particular, can be used to probe experimentally the 
quantum field effect for two-atom cavity quantum electrodynamics. 
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On the other hand, for lower coupling strengths, we can safely ignore 
the effects of rare two-atom events in cavity quantum electrodynamics. 
The theory presented here allows us to distinguish the two cases and 
know when two-atom effects are important. 
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A Eigenvectors and eigenvalues for two 
two— level atoms in as single— mode cav- 



The Hamiltonian for two 2LAs in an optical cavity can be derived 
from equation (J5J) and is given by 



For the one-atom case, the Hamiltonian is diagonalised to yield the 
JC dressed states. Diagonalising equation (|2U)) for g\ = g2 yields the 
dressed states of the two-2LA Tavis-Cummings system |15j . In this 
section, we consider the dressed states for g\ / #2- 

The (one-atom) JC spectrum of figure |2f a) is a ladder of states 
consisting of a ground state and then a sequence of couplets. The 
dressed states are designated as \n) £ , with n designating the couplet 
number and e £ {+, — } designating whether the state is the one with 
higher or lower energy. The quantity n corresponds to the number of 
quanta in the system; for example the states |1)± are superpositions 
of |l)|g) and |0)|e), i.e. a superposition of a photon (and the atom in 
the ground state) with an atomic excitation (and the cavity mode in 
the ground state) . Here we generalise this approach to the case of two 
2LAs in the cavity. 



ity 




(20) 
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The first few levels of the spectral ladder are shown in figure E^b). 
The ground state is designated as |0) and corresponds to an absence 
of photons in the cavity and both atoms in the ground state. For g 2 = 
g\ + g%, the set of one-quantum states form a triplet, and the three 
eigenstates of this triplet are 

|i)o = g-^i&Wls) - gx\s)\*)) (21) 

|1)± = T2- 1/2 i|0)|g)|g)+2- 1 /2 r i| )( 52 |g)|e) + 9l |e)|g)),(22) 

with corresponding eigenvalues 

A ± = to±g (for |1>±) (23) 
A = u (for |1) ), (24) 

respectively. 

The multi-quanta states, for two and more energy quanta, are 
quadruplets. We use the notion |n) ee /, with n the number of quanta 
(n = 2 for the first quadruplet, n = 3 for the second, etc.). The 
subscripts e and e' are each either + or — (i.e. e, e' 6 {+, — }), and 
there are four distinct combinations of e e' corresponding to each of the 
four states in the quadruplet. The choices of e and e' are determined 
by signs in the expression for the eigenvalue of the particular state. 
The four states of the n th quadruplet are 

\n + 1) + ± = Ai n) |n + 1) |g) |g) - i([ n 2\n) |g) |e) 

-iC£ ) |n)|e)|g) + (AAi n) )~ l |n-l)|g)|g) 

|n + 1)_± = Ai n) |n + l)|g)|g) + iCi±V>|g)|e) 

+<£V)|e)|g) + (AAi" ) )-5|n-l)|g)|g). 

The coefficients employ the convenient terms 



:(n) 



(2n + 1)2( 5 2 + g\Y - An{n + l)(gf - g\)\ 



A(B) / _ {gl+gl±^ 



4:gig 2 y/n(n + l) 
{n)' _ ff 2 2 +gf(l + 4n)±»(") 



Ci± 
C 2± 



2 52 ^2n((2n + l)( 5l 2 +5 2 )±H (n) ) 
{n)> _ g 2 +g 2 2 (l + 4n)±»W 



2< 7l ^2n((2n + l)( 5l 2 + 5 2 )±HW) 
17 



with 

^W =A (n)' 2 + c (n)' 2 + c (n)' 2 + 1) 

An)' An)' An)' 

jAn) _ A ± Ap) _ M± An) _ L >2± 

The corresponding eigenvalues of the n quadruplet are 

^ _ („ + !)„ ± ^ + (for |„ + 1) ±± ),(25) 

AJ± = (n + 1)u T y (2„ + i) to? ; g i )± HC) (for |n + i)t±) (26) 

For this analysis, we have assumed that g\ ^ gi- If g\ = g2, we 
obtain the Tavis-Cummings model |15j . for which the ladder consists 
of a ground state singlet, a doublet for the one-quantum states (with 
an energy splitting A2g) and then triplets for two or more quantum 
states. As this case has been studied in depth, we do not analyse the 
Tavis-Cummings model here. 
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Figure 1: The scaled coupling strength distributions nP{g) as a function 
of g/g max for single atoms passing through an optical cavity supporting a 
single TEMqo mode. The solid curve corresponds to a typical distribution 
for a rectangular mask filtering the atomic beam. The dotted line corresponds 
to the absence of a mask. 
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Figure 2: The lowest dressed-state multiplet in the nonlinear regime for 
(a) the Jaynes-Cummings system and (b) two two-level atoms coupled to a 
single quantised field mode of the cavity. Each system may be driven by a 
bichromatic field with frequencies u)\ and u>2- Two excitation pathways to 
the second multiplet are depicted in each case. Inhomogeneous broadening 
of multiplets is depicted by the shaded region. In case (a), the first doublet 
has inhomogeneous broadening 2g max , and the second has width 2\/2g max . 
In case (b), the triplet has an inhomogeneous broadening width of 2\/2g max . 
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Figure 3: Two-photon count rate vs normalised scanning frequency for two 
atoms in a sparse atomic beam for g = gf = 63k, 8\/k = l/y/2, 82/^ = v2 
and j/k = 2 (a) without and (b) with background subtraction. Vertical 
dotted lines are placed in (b) to identify the obvious peaks. 
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Figure 4: The excitation pathway responsible for the two-atom signature in 
2PCS. 
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Figure 5: The 2PCR (a) without and (b) with background subtraction as 
a dotted line. The solid line is the 2PCR for the c^i-driven |0) « — ► |1)_ 
transition artificially suppressed in the simulations. The dotted lines are 
replicates of those in figure a, b), and the scales here correspond to those 
scales. 
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Figure 6: The 2PCR without (first column) and with (second column) back- 
ground subtraction as the dotted line. The solid line corresponds to the 2PCR 
for an ^-driven transition artificially suppressed in the numerical simulation. 
The suppressed transition is |1)_ < — ► |2) ++ for the first row, |1)_ < — ► |2)+_ 

for the second row, |1)_ < > |2) for the third row and |1)_ < ► |2)_+ for 

the fourth row. 
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Figure 7: The 2PCR without (first column) and with (second column) back- 
ground subtraction as the dotted line. The solid line corresponds to the 2PCR 
for an tc^— driven transition artificially suppressed in the numerical simula- 
tion. The suppressed transition is |0) « — ► |1) + for the first row, |0) « — ► |l)o 
for the second row and |0) < ► |1)_ for the third row. 
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Figure 8: The 2PCR without (first column) and with (second column) back- 
ground subtraction as the dotted line. The solid line corresponds to the 2PCR 
for an u;i-driven transition artificially suppressed in the numerical simulation. 

The suppressed transition is 1 1 ) o < — ► |2) for the first row, |l)o « — ► |2)_+ 

for the second row and |1)_ < ► |2)_+ for the third row. 
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Figure 9: The 2PCR without (first column) and with (second column) back- 
ground subtraction as the dotted line. The solid line corresponds to the 
2PCR for an c^i-driven a transition artificially suppressed in the numeri- 
cal simulation. The suppressed transition is |1} + < — ► |2) + _ for the first 

row, |1) + < ► 1 2) for the second row and |1) + « ► |2)_ + for the third 

row. 
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Figure 10: Two-photon count rate vs normalised scanning frequency for one 
atom (dotted line), two atoms (O) and a sparse atomic beam (solid line) 
with P1/P2 = 9 for £\J k = l/y/2, £2/ n = V% and 7//C = 2 with (a) gf = 63k 
and (b) gf — 9k. 
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